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Abstract 

We study the thermodynamical phase structures of holographic QCD with nontrivial topo- 
logically charged domain- wall/membranes which are originally related to the multiple 6'- vacua 
in the large N(, limit. We realize the topologically charged membranes as the holographic D6- 
brane fluxes in the Sakai-Sugimoto model. The D6-brane fluxes couple to the probe D8-D8 via 
Chern-Simon term, and act as the source for the baryonic current density of QCD. We find 
rich phase structures of the dual meson system by varying asymptotic separation of D8 and D8. 
Especially, there can be a thermodynamically favored and stable phase of finite baryonic current 
density. This provides the supporting evidence for the discovery of the topologically charged 
membranes found in the lattice QCD calculations. We also find a crossover phase with the lim- 
iting baryonic current density and temperature which suggest a Hagedorn-like phase transition 
of meson dissociation. 
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1 Introduction 



One of the important issues in QCD is to understand its vacuum structure. The non-perturbative 
vacua such as the ones dominated by the instanton configurations [IJ are relevant in many aspects 
of QCD. For example, it was proposed that the instanton is related to the axial U(l) anomaly and 
yields large r/' mass [2]. However, the effect of dilute instanton gas is e~^'' and negligible in the large 
Nc limit. To account for such a discrepancy, a co-dimensional one object of continuous topological 
charge distribution, namely, the topologically charged membrane was proposed in jSlH], and later 
on realized as the D6-brane in the context of AdS/CFT correspondence [5]. Interestingly, in parallel 
with the Wilson loop summarizing the Schwinger effect in the CP" model, Liischer proposed in [3| 
a Wilson bag term for a 3-form potential C3 over a 3-dimensional hypersurface to characterize the 
pair production rate of the topological domain wall, i.e., 

(e*^^-^^)e=ooce-^^/o"^'^^'^(^') (1) 

where F{9) = —i-^^{tvFF)g is the topological charge density with respect to 0- vacuum, C^,yp = 
—tT{Afj_AiyAp + ^A^^diyAp-\) is an abehan gauge field of fourth kind such that Ae"''^'^^ d^^C /s^yS] = tiFF. 
The 3-dimensional closed hypersurface S with enclosed 4-volume V4 describes the history of the 
pair creation and annihilation of the topologically charged membranes. This is the generalization 
of the Coleman-Schwinger effect for (l+l)-dimensional quantum electrodynamics [6]. 

Typically, a 3-form potential is sourced by some charged membranes and will create a long-range 
Coulomb-type constant force in the ambient (3+l)-dimensional Minkowski space [7]. If the force 
further couples to light particles, then the Schwinger-like pair production will diminish the 3-form 
potential and its production rate can be described by ([1]). However, the 3-form potential C3 is not 
an elementary field but a composite one sourced by the topologically charged membranes, it is not 
clear how it couples to the pair-produced light fields. Later, we will see that the 3-form potential 
will couple to the baryonic current from the holographic Chern-Simon term. 

For finite Nc, the lattice QCD is reliable and practical to explore the non-pertubative QCD 
phenomenon. Implementing the lattice chiral symmetry and the associated topological charge 
density operator, the evidence for the topologically charged membranes in pure-glue SU (3) lattice 
gauge theory has been reported in [8] . Guided by the non-positivity of the two-point function of the 
topological charge density at non-zero distance [9], which rules out the instanton gas, they found 
that the dominated configuration in the vacuum is permeated by two oppositely charged sign- 
coherent connected structures ("sheet"). Each sheet is built from elementary 3D cubes connected 
through 2D surface. 

On the other hand, it is believed that this topologically charged membrane is related to closely 
juxtaposed dipole-like D6-D6-branes in the holographic QCD model [5l \W\ without considering 
its effect on the hadron dynamics. The holographic QCD model has been generalized by Sakai- 
Sugimoto \12\ [T3] to include the hadron physics by using the D4-D8-D8 brane configuration to 

^In the literatures it is usually called topologically charged domain wall. In this paper we prefer to call it 
topologically charged membrane. 
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capture the quark dynamics. This model is based on Witten's model for 4-dimensional pure Yang- 
Mills theory by wrapping Nc D4-branes on a Scherk-Schwarz circle [Tl|, and then putting additional 
probe Nf D8 and D8-branes transverse to the circle. The open strings connecting D4 and D8 or 
D8 then provide the chiral fermions as quarks in the fundamental representation of both the gauge 
group U{Nc) and the flavor group U{Nf) l x U{Nf) /j. In the strongly coupled regime in the large Nc 
limit, the D4-branes are condensed into Witten's geometry [11], which is of cigar-shape so that the 
D8 and D8 are curved and smoothly connected into U-shaped configuration with an asymptotic 
separation L at infinity. The worldvolume theory of the probe D8-brane in Witten's geometry 
successfully realizes the meson and hadron physics [12\ I13j . This U-shaped D8-D8 configuration 
geometrically realizes the spontaneous chiral symmetry breaking at low energy. Besides that, the 
Sakai-Sugimoto model has also been generalized to the deconfined phase [15] . 

Motivated by the appearance of the topologically charged membrane in the lattice simulation [8j 
and its implication of Liischer's Wilson bag, it is interesting to study this effect in the various phases 
of the Sakai-Sugimoto model. In this paper, we will realize the topologically charged membranes by 
the Ramond-Ramond (RR) 8-form fluxes, which can be understood as condensate of the smeared 
D6-D6 branes. Interestingly, from the Chern-Simon term of the probe D8-D8 branes, we find that 
the RR 8-form will couple to the spatial component of the U(l) part of the flavor gauge field, 
which is holographic dual to the coupling between Liischer potential C3 and the baryonic current. 
Moreover, realizing the external baryonic fields as the D6-branes' flux helps us to understand 
its back reaction to the D8-D8 probe branes' shape, and remove the singular chemical potential 
configuration. This makes possible the exploration of the full parameter space of the phase diagram 
of Sakai-Sugimoto model with nontrivial vacuum dominated by topologically charged membranes. 
Similar story has happened in introducing the finite baryon number density in Sakai-Sugimoto 
model |19| . our investigation for relation between the topologically charged membrane and the 
baryonic current will follow the same line. 

Recently, the baryonic current coupled to the external baryonic electric field in Sakai-Sugimoto 
model has been considered in [18], where the baryonic Ohm's law for the decay of the baryonic 
electric field is identified and can be understood as the Schwinger effect. Instead, we find that 
there is no imaginary part in D8's DBI action, and the topologically charged membranes realized 
as the D6-branes' fluxes are dynamically stable. Therefore, the Schwinger effect as suggested in ^ 
or [5] through the nucleation of the baryonic current is absent for topologically charged membrane. 
However, we find very rich phase structures of thermodynamics. 

The paper is organized as following. In the next section we give configuration of the source 
D6-brane fluxes as the holographic dual of the topologically charged membranes in QCD. Especially 
our configuration will couple to boundary baryonic current density and yield non-trivial thermo- 
dynamics. In section 3 the cusp configuration of probe D8-D8 in the confined phase is solved, its 
holographic thermodynamic is derived and the phase structures are displayed. We find that there 
are rich phase structures as we vary the asymptotic separation of D8 and D8. For example, we 
find a crossover phase with a limiting baryonic current density indicating a Hagedorn-like phase 
transition of meson dissociation, i.e. quark deconfinement. More important, there exists stable 
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topologically charged membranes which support the lattice QCD's results in [8j. Similar consid- 
erations for the high temperature phase of holographic QCD are done in section 4 in which there 
are more rich phase structures due to the additional temperature dependence. Finally we conclude 
our paper in section 5. 



2 The topologically charged membranes in Sakai-Sugimoto model 

The holographic dual of 4-dimensional, large Nc QCD with massless flavors was proposed in [121ll3j. 
This model is constructed by embedding Nf D8-D8 probe branes in the D4-brane background. In 
addition, we will also turn on the condensed A'^g D6-branes as the topologically charged membranes. 
The brane setting is arranged as following: 

(4) 5 6 7 8 9 
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The dual QCD lives on i?^ spanned by the coordinates to where the above D6s look like the 
co-dimension one domain wall/membranes. For simplicity, we will only consider the topologically 
charged membrane with x"^ as its transverse direction. Furthermore, we will smear the D6-brane 
fluxes along the x^-direction uniformly, which represents an uniform topologically charged mem- 
brane distribution. In the above brane setting, the D4 and D6-branes are condensed into geometric 
background but with the D8-D8 as the probe brane whose worldvolume theory will capture the 
essence of QCD's meson physics [12]. Moreover, we assume 

A^6 « Nf « iVc, (2) 

to suppress D6s' back-reaction to the D4 background and also its contribution to the vacuum energy 
of the dual meson's (thermo-)dynamics of D8-D8. 

As the supergravity background, the above D6-brane fluxes satisfy the source equatioiJ§ 

duiV^F^'^''-^'^) = J, (3) 

where U is the radial coordinate of Witten's geometry. The explicit form of the D6-fluxes Fg is 
determined by the D6-brane source Jg. To have an energy-scale- independent thermodynamics of 
the holographic QCD, it turns out that we have to introduce the D6-fluxes in the following form 

Fg = Nq 6{U - Uc)dx° A dx^ A dx^ AdU A dVl^ (4) 

where the D6-branes are uniformly smeared along the xs-directiorlf]. We then plug the F^ of dH) 
into source equation ([3]) to deduce the source Jg. From the term di/6{U — Uc) we can see the leading 



^In this paper, we consider field equations in the string frame. 

smooth [/-profile D6-fluxes may introduce inhomogeneous current density profile along [/-direction, this is 
beyond the scope of simple thermodynamics considered in this paper. 
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terms of contain both positive and negative singular charge density at U = Uc- This imphes a 
closely packed D6-D6 pair located at U = Uc- This is consistent with the lattice result in [8] in 
which appears a vacuum with two oppositely charged sign-coherent connected structures ("sheet"). 

Before we continue, we should make sure the D6-fluxes will not cause large back reaction to the 
D4-branes' background metric. To check this, we only need to compare the on-shell actions of the 
D4 and D6-fluxes. The ratio of D6 to D4's kinetic terms (in string frame) is 

^-§m-uc)? (5) 

where the 4- form fluxes F4 is sourced by the Nc D4-branes, and we have omitted the irrelevant 
power of U factor. With the assumption ([2]), the D6s' kinetic term is suppressed if we firstly 

regularize the singular profile dH) and then take large Nc limit. Moreover, in this way the D6-fluxes 
will not add the vacuum energy to the (thermo-) dynamics of the dual meson system of D8-D8. 

Given ^ we can obtain a Hodge dual 2-form F34, which is associated with an axion 04 by 
-F34 := ^304 — 5403. The axion couples to the topological charge density of dual QCD 

04 tiFF. (6) 

At the same time, the 6 parameter of dual QCD vacuum can be defined as 

lim / dx^a4 = / dU dx'^FuA (7) 



Ukk 



which we have used the Stokes theorem. Note that this yields zero value of 9 for our D6-fluxes 
since Fjj^ = by construction. This is in contrast to the case considered in [SJ EH] in which a D6's 
2-form Fjji = n^/U^ corresponding to = in ([3]) is used, so that the 9 parameter is nonzero 
and finite. 

In Sakai-Sugimoto model, the meson dynamics is described by the DBI action of the probe 
D8-D8-branes with U{Nf) gauge fields as the holographical dual of the meson fields. Especially, 
the diagonal U{1) part of the fiavor group is holographically related to the baryonic current ip°'^^ipa 
|181 119j. Besides, we should also include the Chern-Simon term 

Scs = nj Cp+i A tT{2^a'Ff-^^'^ (8) 

if there are non-trivial RR-form fields Cp+i generated by Dp-brane sources. In the presence of 
above D6-branes acting as the topologically charged membranes in the holographic dual QCD, the 
source Chern-Simon term to be incorporated is 

Scs = TsJc7A tr(27ra'F) = iV/Tg J Fs A asdx^ (9) 

where 03 = ^™ is the x^-component of the diagonal U(l) part of the gauge field. As usual, 
can be considered as the chemical potential of the baryonic current. On the other hand, if we 
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identify 04 in ([6]) as the baryo-axio proposed in i.e., F34 = *io-^8) then the suggestive form 
of dS]) provides a holographic reahzation of the conjectured couphng proposed in [7] between the 
baryo-axion and the baryonic current. In some sense, the spontaneous baryonic current induced by 
the topofogicahy charged membranes is analogous to the one by the skyrmions. 



3 Confined phase with topologically charged membranes 

In this section, we will discuss the thermodynamics of Sakai-Sugimoto model in the confined phase 
with topologically charged membranes discussed in the previous section. The main idea of the 
Sakai-Sugimoto model is to realize the meson dynamics of strongly-coupled QCD by probe D8-D8 
branes in the Witten's geometry of Nc non-extremal D4-branes, which is holographically dual to 
strongly coupled large Nc Yang-Mills theory in the confined phase. 

The Witten's geometry of the D4-brane background describing QCD confined phase is given by 

ds' = (I) ' {rj.^dx'^dx'^ + fmdx^) + (I) ' (-^ + U'dnl^ (10) 

where /i, = 0, 1, 2, 3, and Q4 = Svr^/S and 64 are the volume and volume form of the internal unit 
S^. The parameter R is related to the string coupling and string length by R^ = ngsNJ^, and the 
metric describes a throat geometry with the circle smoothly shrink to a tip at the IR endpoint 
U = Ukk by requiring the period of to be 47ri?^/^/3C/j^{^ = 27r/MKK to avoid conical singularity. 
Moreover, the relations between the gravity and the gauge theory quantities are 

where is the gauge coupling of the dual 4-dimensional Yang- Mills theory. 

The Nf probe D8-D8-branes are embedded in such a way with X4-direction as the only transverse 
direction, this yields the induced metric on its worldvolume as follows 



3 , , 3 



dsm = ( Vfifdx^dx'" + 



+ dU' + U'd^l 



(13) 



The D8-D8 worldvolume dynamics is described by the usual DBI action with U{Nf) gauge field 
as follows 



Sdbi = -NfTs J d9xe-<^V-det(5 + 27ra'F) = -N J dUU^^ f{U){dux^y + (^|^ (-^ + (^03) 

(14) 

^Compactifying C7 in ((9)1 on the internal 4-sphere of the Witten's geometry of D4-branes, the reduced 3-forni 
potential can be thought as the 4-dimensional Hodge dual of the baryo-axion, that is the field sourced by the 
topologically charged membrane. 



2^ 



5 




where we have defined N = ^/'^s^'tVa/j ^j^g^g V3 = / dx^dx'^dx^ and /? = / dx^, and the x^- 
component of the diagonal U(l) part of the gauge field 03 = ^'-^ . 

Prom the dictionary between the bulk and boundary quantities, for example the ones given in 
[19] for the discussions of the finite holographic baryon number density, the boundary value of the 
above bulk background gauge field 03 corresponds to the chemical potential of the x^-component 
of the baryonic current of the dual QCD. The dual baryonic current density can also be identified 
from the normalizable mode of 03. 

Moreover, we should also add the Chern-Simon term ([8]) for the nonzero RR 8-from ([4]) corre- 
sponding to the topologically charged membranes in the QCD vacuum. This yields 

Scs = ^nhj dU6{U - UMU) (15) 

which acts as a source term for the boundary baryonic current. In the above, we cut off the volume 
of the x^-direction to L3, and define the baryonic current density of the dual QCD by 



NjTsN.L, V,= fdx^dx^ (16) 

N 87r2 V2f5 ' ^2- J ax ax . Ubj 



Solving D8-D8 configuration with a cusp 

The Euclidean on-shell action of the probe D8-brane will encode the thermodynamics of the holo- 
graphic QCD. To acquire the information, we need to solve the D8-brane configuration. 

In flat space-time, the D8-brane will be straight and flat. However, due to the cigar shape of 
Witten's geometry, the probe D8-brane will force to bend into D8-brane at Uq. Integrating ([20]l to 
get x^{U) with fixed asymptotic separation of D8 and D8, one will get a U-shaped D8-D8-brane 
configuration, which presents the chiral symmetry breaking of the dual QCD in a geometric way 
|12j . Moreover, once the D6-brane fluxes are turned as ([Ij), the above smooth U-shaped D8-D8 will 
then develop a cusp as shown in Fig(T] due to the pulling of the force from the D6-brane source, 
which signifies the existence of a new thermodynamical phased. 



^This is similar to the case of finite baryon number density, for example see the last reference in [19] . 
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We first solve the gauge field 03 from its equation of motion derived from the total action of 
D8-D8, which is 

= N Q{U - Uc), Uc > Ukk (17) 
where the properly normalized conjugate momentum to 03 is defined by 

_ SSpBi _ UR^dyas 

'-'-as = TT^ T - I o UoJ 

^^^'''''^ ^fmdux,r + {§f{j^ + {duasr) 

and the step function on the right-handed-side comes from the delta function D6-fluxes source 
encoded in (fT5]) . 

From ([18]) we can get 



UR3' 



. „g \ (19) 

The above choice of the overall sign is to be consistent with (|17p and p8|) . 

On the other hand, there is no source term for the equation of motion of x'^, we instead need 
introduce an integration constant C. With the help of (jl9p . we have 

with 

HiU) ^ U^fiU) (1 - ^) . (21) 
The integration constant is assigned by relating it to the proper cusp angle 



tan^c = .l-^dux'^ \u=Uc (22) 

guu 



and it yields 



C^C/cm)(l-^)sin2 0c. (23) 



To further integrate out (I20p to obtain the full D8-D8 configuration we need one more condition, 
which is to fix the asymptotic separation between D8 and D8 

L = 2 dU dux^. (24) 

JUc 

Moreover, the size of the proper cusp angle 9c can be determined by the force balance condition 
between the D8-D8's tension /ds and the pulling of D6 /dgj i-e., 

/d8 + /d6 = 0. (25) 

The proper tension of D8-D8 is obtained by varying its Hamiltonian density with respect to the 
proper distance Uc along the D8 brane, that is, 

1 5Hds I 



1 



nrLyUc) \ I r-,^ \U=Ur, with ni,,L fixed 

Jjj^ 6{dux^) dU 



4^ 

(27) 
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where the Hamiltonian density is the Legendre transformation of the DBI action ()14p as 



Hd& = -Sdbi + dU Ua.duas 



with 



n {x\u)) = 



\ 



f{U)idux^)^ + 



nt 



f{U) 



(28) 



(29) 



The resultant proper D8-D8 tension is 



i 



riT 



cos ( 



(31) 



To obtain the pulhng force of the D6-fluxes, we need to know the on-shell Chern-Simon action 
which can be obtained by plugging (f20l) into (fT9]) . and then integrating it over U to get 03 (C/). 
Result antly, 



Scs\on-shell = ^^^603(^^0) 

Therefore, the force by D6 flux is 
/d6 = 



dU- 



n 



my 



Uc 



VH{U) - C 



(32) 



1 5Scs 



nt 



I ;tt \U=Uc, with nh,L fixed~ p \ n 

Given the explicit form of Jds and foe, the force balance condition ([25|) yields 



(33) 



COS^ 9r 



WW 



ufW 



(34) 



Absence of dynamical instability via tunneling 



Though the factors H{U) — C and ^ 1 — jj^w a-PP^aring in the on-shell action seem implying 
the possibility of dynamical instability via tunneling such as the Schwinger effect if these factors are 
turned into complex. This usually happens for a DBI action with the electric field larger than its 
maximal field strength. However, recall that we in fact turn on not electric but magnetic field dua^ 
for the holographic baryonic current density. Thus, we will not expect such a dynamical instability 
in our case. This is also in contrast to the case considered in jT9] for the finite holographic baryon 
number density, in which a electric field Fuo is turned on in the DBI action, so there is a possible 



While varying dux"^ with Uc we should keep the asymptotic separation (|24|l fixed, this yields 



dux\Uc) = 



Uc 



djdux") 
dU 



\U — U ciWith nf,,L fixed 



dU 



(30) 



from which we can evaluate the second term in (|27p along with the Hamilton's equation du staj^x'^) ) ~ ^' 
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dynamical instability for the Lorentzian configurations. Indeed, (j34p implies ^75^ < 1/2 so that 
cos^ 9c < 1 for a physical cusp configuration, thus the above factors are real for U > Uc and 
the dynamical instability is absent. This implies a maximal baryonic current density rif, < a / — ^ — 
beyond which the cusp configuration cannot be sustained, and it could imply the meson dissociation. 

Thermodynamics of holographic QCD 

From the usual thermodynamical and holographical interpretation of the Euclidean path integral 
in the grand canonical ensemble at temperature T 

n 

we can identify the on-shell Euclidean total action as the Gibbs free energy density with fixed 
chemical potential i.e., 

on— shell ~l~ ^Cslon— shell- 

(36) 

The holographic chemical potential of the baryonic current in dual QCD is identified as as the 
boundary value of gauge field, that is, 

/i = No3(oo) (37) 

and the on-shell Euclidean Chern-Simon action differs from (j32p by an overall minus sign, together 
with the on-shell Euclidean DBI action obtained from (1141). (1191) and (1201) 



SDBl\on-shell = N / dU —===. (38) 

JUc ^/H[U)-C 

Instead, we can also work in the canonical ensemble with fixed charge density by performing 
the Legendre transformation to obtain the Helmholtz free energy density 

F{T, m) = E -TS = n + fini, (39) 

and from (j32p and (j38p we obtain 

3 13 / \ 



which is nothing but the on-shell value of the Hamiltonian density (j28p and (129p as expected from 
the definition of the Helmholtz free energy 0. Moreover, the relation between chemical potential 
and charge density in the canonical ensemble can be obtained vij§ 

Km) = (41) 



^In the confined phase, the temperature dependence is trivial and its thermodynamics will be thought as the zero 
temperature case. 

*Or, we can vary the Hamiltonian density (|28p so that 

, s SHds I dUc , dHns , 777T, , (J-ffog 
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Inverting this relation we can get f^(/u) = i7(/U, Furthermore, the thermodynamical stabihty 

is characterized by the baryonic current susceptibihty 

If X is positive, it is stable; otherwise, it is not. 
Phase structures 

Since the condition (j24p of fixed asymptotic separation is nontrivial so that it is impossible to 
obtain the analytical form for the above thermodynamical quantities, instead we will substitute the 
proper cusp angle by (j34p and numerically solve Uc{nt,) from (j24p and obtain the thermodynamical 
quantities. Moreover, we may expect rich phase structures due to the factors such as y/ H{U) — C 
in the thermodynamical quantities. This is in contrast to the simple phase structure considered in 
|19j for the case with non-zero holographic baryon number density. In fact, the most peculiar we 
have found is that the thermodynamical behaviors strongly depend on the value of the asymptotic 
separation L. As L is varied, there are the following phases: (a) thermodynamically unfavored and 
unstable phase, (b) thermodynamically favored and stable one, (c) thermodynamically unfavored 
but stable one, and (d) crossover ones between the above phases. The feature of the rich phase 
structures is quite similar to the one of holographic entanglement entropy characterized by a critical 
asymptotic separation [T6j. 

In our numerical computations we have set [§ 




i?3 

Ukk = 0.5, R = 1^L<—\I- — - 2.96. (43) 

Ukk 

Moreover, the thermodynamics for the confined phase for a given L can be just summarized in 
two plots, one is ^ v.s. rib, the other is Afi v.s. n, from which we can obtain AF{nb). The AO 
is the difference of the Gibbs free energy between topologically charged membrane state (rib 7^ 0) 
and vacuum state {nb = 0). Similarly for AF(nfe). It is also interesting to see how the tip of cusp 
changes while varying Ub from Uc{nb), however, we omit it in the following because it is not relevant 
for thermodynamics. 

Thermodynamically unfavored and unstable phase 

The results of a typical case with L = 0.4 for such a phase are summarized in Figl2][3j We can see 
that: 



and obtain ^f- by varying the condition (I24p over nt so that 



This is useful for numerical implementation. 

^We should caution the readers here. In this chosen set of units, the 't Hooft coupling is not much larger than one 
so that the SUGRA approximation is not good. However, it is conventionally used for the numerical computations 
in the literatures [181 119] with surprisingly good results in comparison with experimental value. 
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• AQ in Figl3]is positive so that the phase is thermodynamically unfavored. 

• Moreover, this phase is also thermodynamically unstable as can be seen from the negative 
susceptibility derived from the slope of FiglJl This is also suggested by the negativity of the 
chemical potential fi. 

• In these plots, only the part with the baryonic current density rib > 0.2 is shown. In our 
numerical results not shown in Figl2]l3l we find that for very small the Gibbs free energy is 
far higher and is thermodynamically unfavored. Therefore we just show the parameter space 
with smaller Ail. region which is also more realizable without numerical convergence problem. 
Similar situation happens for the other cases of different L's. 

Crossover phase 

In between the thermodynamically unstable phase for L < 0.6 and the stable ones for L > 0.7, 
there is a crossover as shown in Figl3]l5]for L = 0.65. We see that 

• Two branches appear in FigU] because there are two values of chemical potential for a fixed 
baryonic current density. 

• One branch is still thermodynamically unfavored and unstable, but the other is both ther- 
modynamically favored and stable as can be seen from FigU] and El 

• From FigHlwe see there exists a maximal value of Ub at which the both chemical potential and 
susceptibility are divergent, and thus limiting. This is similar to the limiting susceptibility 
for the Hagedorn phase transition [20\ I14j. However, this limiting baryonic current density 
is smaller than the one required by cos^ 6c < 1 to sustain the cusp configuration. This 
suggests that the cusp configuration cannot be sustained for the smaller rih than expected, 
and the limiting behavior might be a signature for meson dissociation which bears Hagedorn- 
like behavior as expected. This could be related to the similar Hagedorn behavior for the 
holographic entanglement entropy considered in |16j . 




An 
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-I.I 



-1.0 



Figure 2: v.s. Ub {L = 0.4) 



Figure 3: AO. v.s. ^^ {L = 0.4) 
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Figure 4: ^ v.s. Ub (L = 0.65) Figure 5: AO v.s. n {L = 0.65) 



Therefore, our result shows that the topologicahy charged membranes can induce quark de- 
confinement even at zero temperature. 

• Another interesting feature for this phase is the sign change of the chemical potential for 
baryonic current density as we vary L. This is different from the usual crossover in QCD 
of finite baryon number density, for which there is no sign change of chemical potential for 
baryon number density [21j. Instead, it reminds us the BCS-BEC crossover in condensed 
matter system [22j where there is a sign change of chemical potential from the BCS phase 
(// > 0, weakly coupled) to the EEC phase {fj, < 0, strongly coupled) as the number density 
increases. 

Thermodynamically favored and stable phase 

One of the goal of this work is to search for the new stable vacuum of large Nc QCD with non- 
trivial topologically charged membrane structure. The results of a typical case with L = 1 for such 
a phase are shown in Figl6][71 

• Ail. in Figl7]is negative for the most of positive fi region, it implies this phase is thermody- 
namically favored. It is interesting to see there is some unfavored region for negative as the 
tail of the crossover. 

• Moreover, this phase is also thermodynamically stable as can been from the positive sus- 
ceptibility derived from the slope of Figl6l Similarly, there is a small region of negative 
susceptibility as the tail of the crossover. 

• Therefore, there exists a thermodynamically favored and stable phase with non-zero baryonic 
current density in the confined phase of holographic QCD. This provides the supporting 
evidence to the finding of the topologically charged membranes in the lattice QCD calculations 
reported in [8]. 
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Thermodynamically unfavored but stable phase 

As L goes higher, the positive AO part will gradually dominate the negative one, and finally we will 
end up at such a phase where the negative AO part will disappear but it is still thermodynamically 
stable. The typical results with L = 1.5 are given in FigHJlHl and this phase persists for larger L 
0. It is straightforward to read off the thermodynamical behavior as before. Since AO is small, 
it suggests that the topologically charged membranes may have quite a chance to sustain in the 
vacuum even it is relatively unfavored. 

In summary, as the asymptotic separation L of D8-D8 increases, the meson system of the holo- 
graphic QCD goes from a thermodynamically unfavored and unstable phase through a crossover, 
and then becomes thermodynamically favored and stable. Finally, it ends at a thermodynamically 
unfavored but stable phase. Such a rich phase structure at zero-temperature system is expected 
for a strongly interacting theory as holographic QCD. This is in contrast to the free theory with 
trivial thermodynamical phase at zero temperature. 

^"Our numerical calculation is reliable up to L « 2.2, for higher L we have serious convergence problem since Uc 
is very close to ?7kk and it becomes harder to solve p4p numerically. 
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4 High temperature phase with topologically charged membranes 



The background geometry dual to high temperature phase of QCD is obtained from the metric 
(jlOP by wrapping the Euchdean time and cycles as in the usual Witten's cigar like geometry for 
BTZ black hole |1H 114) . This was considered first in [15] and the resultant (Euclidean) metric is 



where f{U) = 1 



3 

I j ' {fiU)dT^ + 5ijdx'dx^ + dxl) + 



Uj \fiU) 



(44) 



[/3 

j^, and the period of the Euclidean time direction is determined to be 



4tt E? 




1 

T 



(45) 



which is identified as the inverse temperature of the dual QCD. As for the confined phase, the 
induced metric on the probe D8-D8 in the background (jl^ is 



D8 



, 3 



- ) {f{U)dT^ + 5i,dx'dx^) + ( - 



, 3 



u 



1 



dU^ + u^dnj 



(46) 



It then results in the DBI action R of D8-D8 



0, 



Sdbi = -^ I dUU^Jf{U){dux^Y+l^^^ {l + fiU)iduasy)) 



(47) 



where N and are defined as for the confined phase. However, the Chern-Simon action is the same 
as p^. 

Since the [/-direction is now no longer caped as in the confined phase but the r-direction is, 
it is then possible for the D8-D8 to disjoin as temperature goes high enough. This indicates the 
chiral symmetry will restore in the high temperature [15]. With this mind we solve the equations 
of motion. 

The equation of motion for 03 is 



dua3 



rib 



and the one for is 



with 



{dux 



4\2 



G{U) - Ct 



G{U) ^ [/« /([/) 



nr 



(48) 

(49) 
(50) 



^^Dilaton's and 4-form flux's profiles are the same as in the confined phase. 

^■^The Lorentzian DBI and Chern-Simon actions are considered until we consider thermodynamics. 
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and the integration constant Ct is related to the proper cusp angle 9c as follows 

CT^Ul{j{Uc)--^^sh^^9c. (51) 

For the chiral symmetry breaking phase the proper cusp angle 9c is nonzero and will be determined 
as in the confined phase. For the chiral symmetry restoration phase, = so that dux'^ = 
Ct = which describes the parallel D8 and D8-branes. Formally, we can consider both cases in a 
unified way, and should be carefully choose 9c and the corresponding quantities when we tune the 
temperature or Ut in numerical calculations. 

The Hamiltonian density derived from the DBI action by the Legendre transformation is 



Ho, = » [iU U']^ {i8,^f + (i)' 7m) - <^^' 

from which we can obtain the proper D8-D8 tension for the cusp development, that is 



3 13 / ffi 

fD8 = -Ni24 Uc' JfiUc)- jj^ cos 9c. (53) 

Similarly we can derive the pulling of D6-fluxes from the Chern-Simon action, and arrive the 
following force balance condition for the cusp (in the chiral symmetry breaking phase) which yields 

cos' 9c = ■ (54) 

Along with the condition for the fixed asymptotic separation of D8 and D8, i.e. (j24p . we can solve 
Uc{nb) as in the confined phase. Note that requiring cos^ 9c < 1, the current density is bounded as 

nb < Uc R^sJ — for a given T, or a bounded temperature as T < -^{Uc — jjt)^ for a given 

As for the consideration of the thermodynamics of the dual QCD in the grand canonical en- 
semble, the Gibbs free energy is the on-shell Euclidean total action and we obtain 

^}{T,n;nb) = -^inb + F{T,nb) (55) 

where the chemical potential fi = No3(oo) as before, and the free energy density 

.oo Rluf(f{U)-j^) 
F{T,nb) = N dU ^ ' (56) 

^f{u) {G{U) - C|) 

is the same as the on-shell value of Hamiltonian density (j52p as expected; however, the nontrivial 
temperature dependence is implicit in the factor f{U) = 1 — which is related to T via (|45p . 
Also, one can obtain /i(T, nb) through 

dF{T,nb) 

n{T,nb) = . (57 

aub 
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Figure 12: An^y v.s. /i^ = 0.7, T ~ 0.207) Figure 13: Afi^p v.s. /i^ = 0.7, T ~ 0.207) 



Phase structures 

At high temperature there are three phases: cusp, parallel (chiral symmetry restored phase), and 
vacuum (rif, = 0). To determine which phase is preferred, we have to compare their Gibbs free 
energies. We denote the Gibbs free energy for cusp phase as f^c, for parallel one as fip and for 
vacuum as ily, and their differences as Aflcv = — ^v, Ailcp = — and A^lpy = 0,p — Qy, 
respectively. Also, the chemical potentials for the cusp and parallel phase are denoted as /ic and 
Hp respectively. 

In order to perform the numerical computation, it is easier to fix either the temperature T or 
the current density as what we will do in the following. We will set i? = 1 as in the confined 
phase. 



Fixed T 

With Ut fixed, the equations of motion are formally similar to the confined phase, we may expect 
similar thermodynamical properties as for the confined phase. Indeed, this is the case and there 
exists a thermodynamically favored and stable phase in the configuration space. Therefore, for 
conciseness, we will not show all the plots but only the relevant ones, and summarize the main 
features as following: 

• For small L (< 0.7 with T ~ 0.169), we find that AQcv, Ailcp, AQpy are all positive, and 
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/Xc, fJ'p are negative. Thus, the vacuum phase is the most favored. Also from the sign 



dlic{n,„T) dfj,p{nt,T) 

thermodynamically unstable. 



of the susceptibility derived from and g^^' , the cusp and parallel phases are 



• As L is rising, there is a crossover (starting around L = 0.7 with T ~ 0.169) and a new 
branch for both AQcv and A^lcp appears as in the confined phase but not for Aflp^. The new 
branch has negative AQ^v and AQ^p so it suggests that the cusp phase is thermodynamically 
favored. Moreover, is positive for the new branch so that the cusp phase is also 
thermodynamically stable. On the other hand, there is no corresponding new branch for 
'^'^^dnl'^^ implies a thermodynamical instability of the parallel configuration. 

• As L goes higher, we find that there appears no new phase beyond crossover, instead there 
exists a critical value of L beyond which there is no physical cusp configuration for finite 
TT-fc. This critical value is T-dependent, for example, ~ 1.28 for T ~ 0.131, L^. ~ 
0.99 for T ~ 0.169, and Lc ^ 0.82 for T ~ 0.207. Besides, the new branch will tend to 
be thermodynamically unfavored and the chemical potential of new branch will tend to be 
negative, although the susceptibility is still positive. This feature is very different from the 
confined phase where the crossover will end with a dominant stable phase after some value 
of L. 

• For higher T, for example L = 0.7 with T ~ 0.207, we find the similar crossover as for 
T ~ 0.169 case, however, we find that A^lcv turns to be positive so that the vacuum phase is 
now the most favored. For more detailed feature, see Fig |10lfT3l Also, as can be seen above, 
the critical value decreases as T increases. 

• Similar to the confined phase, (though not shown here) there is also a Hagedorn-like phase 
transition during crossover with a limiting baryonic current density. 
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Figure 18: n^v v.s. T [L = 0.7, Ub = 0.09) Figure 19: Hc v.s. T {L = 0.7, = 0.09) 



Fixed rih 

It is interesting to observe the temperature dependence for the fixed baryonic current density. The 
main features we will summarize in the fohowing: 

• For small Uf, (e.g. ni, = 0.005), the cusp phase is thermodynamically unfavored. However, 
Ai^cv is quite small so that the topologically charged membrane has quite a chance to sustain 
itself in the vacuum (See Fig |14lfT5]) . So is the parallel configuration. 

It is also interesting to see in Fig|16l there is a critical temperature above which the parallel 
phase can dominate over the vacuum one, and the chiral symmetry is restored. 

• The thermodynamic relation — y = {^)v,u fixed relates the chemical potential to the entropy 
change, where U, V are the internal energy and volume respectively. Therefore, negative 
yields positive entropy changing rate, and implies the instability. 

The chemical potential fic in Figl 171 reaches its local minimum at some temperature beyond 
which it grows but never becomes positive for a physical cusp configuration. That means, 
the above local minimum of fic yields largest entropy changing rate, and relatively unstable 
at some specific temperature. On the other hand, when rn, > 0.08 (see Fig J19p . the chemical 
potential is either monotonically decreasing (old branch) or increasing (new branch), so there 
is no such a special temperature. 
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• For Hb larger than 0.08, a crossover phase appears as shown in Fig llStfTOl with a new thermo- 
dynamically favored and stable cusp branch as seen from AQcu < 0, AQcp < and Hc > 0. 
This new branch is also found in the the fixed T case. 

Moreover, in FigllQI there is a limiting temperature where ^ diverges. This is again of 
Hagedorn-like behavior and implies the meson dissociation as discussed before. This decon- 
finement temperature found here is very close to the chiral symmetry restoration mentioned 
above Ffl. 



5 Conclusion 

In this paper we investigate the QCD vacuum with non-trivial topological charge distribution from 
its holographic dual theory. We find that there exists thermodynamically favored and stable phase 
of mesons in the presence of the topologically charged membranes in both low and high temperature 
regimes. This provides the supporting evidence for the observation of the same topological charge 
structure in the lattice QCD calculation of [S]. we hope our work will inspire more searches for 
this new vacuum in the lattice calculation, and further estimation of its effect on the quark-gluon 
plasma created in the heavy-ion collider. Moreover, from the Chern-Simon term of the probe D8- 
D8-branes, we find that the topologically charged membranes play the role of baryo-axion, which 
couples to the baryonic current density. This is a kind of holographic realization of the baryo-axion 
proposed in [7] . We have also shown that the Lorentzian configuration of our topologically charged 
membrane is dynamically stable, it is interesting to study the modification of the meson and baryon 
spectra in this new QCD vacuum. 

We should emphasize that our D6-brane fluxes used as the holographic topologically charged 
membranes are quite different from the ones for 0- vacuum proposed in [5]. Our D6-brane fluxes 
are sourced by the opposite charged sign-coherent sheet, however, the ones in are not. It is 
interesting to point out that the D6-brane fluxes of [5] will couple to Ajj via Chern-Simon term 
of D8-D8 branes. The lowest KK mode of the gauge field Ajj is identified as pion [12], however, 
this source term will not induce any homogeneous chiral condensate since there is no -F(/4 term in 
the DBI action of D8-D8, which is the same reason for the masslessness of the pion. Instead, the 
Fui^'t' = 1)2,3 in the DBI action of D8-D8 along with the Chern-Simon source term could induce 
inhomogeneous chiral condensate characterized by Au{x^) and the non-trivial 0- vacuum. Further 
study is needed to make sure if the configuration is thermodynamically favored or not. This is in 
contrast to the recent proposal for the mass deformation of the Sakai-Sugimoto model 

Besides, our results show that the thermodynamics of the topologically charged membranes 
strongly depends on the value of asymptotic separation of D8-D8, and a crossover happens around 
some critical asymptotic separations. This critical behavior is quite similar to the one for the 
holographic entanglement where the authors in [16] argue this is related to the Hagedorn spectrum 
of non- interacting bound states. This is in accordance with our interpretation of the Hagedorn- 

^^The deconfinement phase of finite temperature Sakai-Sugimoto model was first discussed in [15j as L is varied, 
however, there is no Hagedorn-hke behavior shown as done here. 
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like meson dissociation phase transition characterized by the hmiting baryonic current density or 
temperature. Despite that, we are not sure why such a critical asymptotic separation exists in 
our case, it may deserve further study. Also it is important to understand the meaning of the 
asymptotic separation in the field theory side. 

Finally, we find there is a sign change of chemical potential as we vary the asymptotic separation 
of D8-D8 around the crossover which also happens in the BCS-BEC crossover of some condensed 
matter system. This analogue motivates the further study to understand the microscopic relation 
between the topologically charged membrane and the baryonic current density and see if there are 
corresponding BCS and EEC phases. 
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